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We have extended the calculation of the correlation functions of heavy quarkonium hybrid op- 
erators with various J quantum numbers to include QCD condensates up to dimension six. In 
contrast to previous analyses which were unable to optimize the QCD sum-rules for certain J^ , 



recent work has shown that inclusion of dimension six condensates stabilizes the hybrid sum rules 

and permits reliable mass predictions. In this work we have investigate the effects of the dimension 

^ . six condensates on remaining channels. After performing the QCD sum rule analysis, we update 

'^'h ' the mass spectra of charmonium and bottomonium hybrids with exotic and non-exotic quantum 

VQ \ numbers. We identify that the negative-parity states with J^ = (0, 1, 2) '", 1 form the lightest 

hybrid supermultiplet while the positive-parity states with J = (0, 1)^ , (0, 1, 2)^"^ belong to a 

heavier hybrid supermultiplet, confirming the supermultiplet structure found in other approaches. 

The hybrid with J^'^ = 0^~ has a much higher mass which may suggest a different excitation of 

T^ . the gluonic field compared to other channels. In agreement with previous results, we find that the 

C*H' J^ = V^^ charmonium hybrid is substantially heavier than the X(3872), which seems to preclude 

f-« ' a pure charmonium hybrid interpretation for this state. 
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I. INTRODUCTION 

In the constituent quark model, hadrons are described as qq mesons and qqq baryons. Most of the expcrimentahy 
observed resonances can be accommodated in the quark model fl|, [2|. However, QCD itself allows a much richer 
hadron spectrum, such as multiquarks, hybrids, glueballs, etc... [3|. Hybrid mesons (qgq) are composed of a color- 
octet quark-antiquark pair and an excited gluonic field. Since the excited gluonic field could carry quantum numbers 
other than 0^^, hybrids can appear with both the ordinary (the same as conventional qq mesons) and exotic quantum 
numbers {J^^' — , 0"' , 1 ^, 2^ , . . .). The states with the exotic quantum numbers are not accessible for a qq 
state. The observation of exotic hybrid states is considered to be an unambiguous signal of gluonic degrees of freedom 
in the hadron spectrum. 

The light hybrids were studied in the MIT bag model [3, Q, in which the hybrids with J^'^ = (0, 1, 2) ^, 1 were 
predicted to form the lightest hybrid supermultiplet consisting of a S-wave color-octet quark-antiquark pair coupled 
to an excited gluonic field with Jg ' ' = 1^ . A higher hybrid supermultiplet composed of a P-wave qq pair and the 
same gluonic excitation would contain states with J^*^ = '", (1 ^)^, (2 ^)^, 3 '", (0, 1, 2)++ where the superscript 
denotes the number of such states 0,01 • Many other methods such as the flux tube model |8l-[l2l|. lattice QCD [l3l - [l6J 
and QCD sum rules |17H24j were also used to study the light hybrids. To date, there is some evidence of the exotic 
light hybrid with J^c = 1-+ [llii. 

For heavy quarkonium hybrids, which hereafter we shall refer to as heavy quar k hybrids, calculations including 
the constituent gluon model [231, the flux tube model [l^l, QCD sum rules jSOl - ISq and lattice QCD jg. [37H41| have 
been performed. The hybrid supermultiplet structures described above were confirmed for the heavy quark sector in 
lattice QCD [y| and the P-wave quasigluon approach [42|- So far, however, no definitive experimental signal for heavy 
hybrid mesons has been observed, although many unexpected charmonium-like and bottomonium-like states have 
been discovered in the past several years p3l - |47| . Some of these states, called X, F, Z, do not fit in the conventional 
quark model easily and are considered to be candidates for exotic states beyond the quark model, such as molecular 
states, tetraquarks, bar yon ium and quarkonium hybrids. For example, K(4260) was interpreted as a charmonium 
hybrid state in Refs. [3, H^ H^]. In Ref. [501, X(3872) was also proposed to be a charmonium hybrid. To search for 
hybrids in the heavy quarkonium region further theoretical investigations of the hybrid spectrum are still needed. 

QCD sum rules is a very powerful non-perturbative method |5ll453l | which has been widely used to study hadron 
structures. In Refs. |30l432| . Govaerts et al. performed a QCD sum- rule analysis of the heavy quark hybrids including 
leading-order contributions up to the dimension four gluonic condensate. The sum rules of some J^'~' channels were 
unstable and thus these mass predictions were unreliable. Recently, the J^*^ = 1~~ [3J], 1"*""*" [3a] and '" [3g| 
channels have been re-analyzed by including the dimension six tri-gluon condensate. The dimension six contributions 
have proven to be very important because they stabilize the hybrid sum rules and allow mass prediction to be made. 
In this work, we will study the mass spectrum of heavy quarkonium hybrids using the QCD sum rule method. We 
will calculate the dimension six condensates including (g^fGGG) and the light-quark condensate (g^jj), where the 
latter condensate has not previously been calculated in the hybrid sum rules. 

The paper is organized as follows. In Sec. II, we calculate the correlation functions and spectral densities of the 
hybrid operators with various quantum numbers and collect them in Appendix|Xl In Sec. HI, we perform the numerical 
analysis and extract masses of charmonium and bottomonium hybrids. In the last section we summarize our results 
and comment on their implications for heavy quarkonium spectroscopy. 



II. LAPLACE SUM RULES FOR THE HEAVY QUARKONIUM HYBRIDS 

We consider the two-point correlation function 

n^.(<z) = * fd^xe^^-^ {0\T[J^{x)4m\0), (1) 



where J^ is the interpolating current. In this work, we use the following hybrid interpolating currents coupling to 
various quantum numbers: 

J^. = 3.0Y7''75G^,g, J^^ - l+-,0— , (2) 

V = 5sQ Y<75GLg, J^^ = 2-+, 1++, 1-+, 0-+ , 



in which Q represents a heavy quark (c or &), gs is the strong coupUng, A° are the Gell-Mann matrices and G° 

is the gluon field strength. By replacing GJJj, with G"^ = ^e^,ya/3G"^'°, we can also obtain the corresponding 
operators with opposite parity. These hybrid correlation functions were originally studied in Refs. |30l432| in which 
the perturbative and gluon condensate contributions have been calculated to perform the QCD sum rule analysis. 
According to these calculations, the stability of the sum rules depends on the relative sign of the gluon condensate 
and the perturbative term. Only a relative negative sign results in stable sum rules. As shown in Refs. J30M32| . the 
sum rules for J^^ = — , 0++, 1+", 1++, 2++ channels are stable while the J^^ = 0"+, 0+^,1"+, 1 — ,2"+ channels 
are unstable. In this work, we will calculate the dimension six condensates in the correlation functions and show 
that they stabilize the hybrid sum rules for all unstable J^'~^ = '",0^' ,1 '", 1 ,2 ^ channels, confirming the 
stabilizing effect of dimension six condensates previously found for 1 [34[, 1++ [3a] and ^ [36| . 
For the hybrid operators in Eq. ([2]), the two-point correlation functions have the following structures: 



I d''xe''^--{0\T[J^{x)4m\0) = 



2 SfJ.l' 



Iiv{q^) + ^-i^Ils{q% (3) 



I j d^xe'^-- {Q\T[J,,{x)Jl{Q)m -\ 



T]fj,pr]iya 'T t]lia1li^p ^Vl-ii^Vpcr 



IlT{q'-) + ..., (4) 



where ry^ii/ = q^q^/q^ ~ g^i^- ny(g^), 115(5^) in Eq. (J3|) and liTiq^) in Eq. (JH) are the invariant structures referring to 
pure spin-1, spin-0 and spin-2 states, respectively. In Eq. (|3|), the invariant structures for spin-0 and spin-1 are not 
written out explicitly because we will not consider contributions arising from these terms in this paper. 
At the hadron level, the correlation function can be described by the dispersion relation. 



^i<l') = i'iT r ,.(/^g ds + '^fKiqr. (5) 

J4m2 s (s - q -te) j^ 



The summation on the right hand side represents the subtraction terms which can be removed by taking the Borel 
transform of Il{q^). The spectral function p{s) is defined using the pole plus continuum approximation 

p{s) =Y,S{s - mi){0\Jp\n){n\jl\0) ^ fxm%Sis - ml) + contimmm, (6) 

n 

where the intermediate states n must have the same quantum numbers as the interpolating currents J^, mx denotes 
the mass of the lowest lying resonance and the dimensionless quantity fx is the coupling of the resonance to the 
current 

{0\Jp\X)^fxmj,e^, (7) 

(0|V|X)-/A-r746,,,. (8) 

At the quark-gluon level, the correlation function can be computed via the operator product expansion (OPE). For 
heavy quark hybrids, the quark condensates are expressed in terms of the gluon condensate via the heavy quark mass 
expansion and hence give no contributions to the correlation function. To calculate the correlation function, only the 
perturbative diagram (Fig.[T|), gluon condensate diagram (Fig.[2|) and diagrams for dimension six condensates (Fig. [3]) 
are involved. In Fig[3l the dimension six condensates contain (g^fGGG) (Fig. [Hb) and (DDG) (Fig. [3K). Using the 
equation of motion, Bianchi identities and commutation relations, the condensate (DDG) can be expressed in terms 
oi {g^JGGG) and {gtn). 

We have calculated the Wilson coefficients using two independent techniques, finding complete agreement between 
them. First, we have used the same approach that was utilized in Refs. 
respectively. In this approach, we utilize the Mathematica package Tarcer 



35, [3g| for '" and 1++ heavy hybrids, 
5JI, which implements the d dimensional 
recurrence relations developed in Refs. [55|,|56[ to reduce the number of distinct integrals that must be calculated. The 
resulting minimal set of loop integrals can be evaluated using results given in Refs. [57 ^], leading to very compact 
expressions for the Wilson coefficients in terms of generalized hypergeometric functions. Second, we have utilized the 
technique used in Refs. |60l - l63j . With this method the perturbative spectral densities are expressed as integrations 
that can be performed numerically. In this approach we use the momentum space quark propagator 



i6ab ^ i_ A^^„ a^"(|^ + m) + (/ + ^W^- 
i) — m 4 2 '"^ (p^ — m^y^ 



iSab{p) = -^^ + -gs^Gl, '' , ' \^^ , (9) 



where a^^ = | [7^ , 7"^], and a , b arc color indices. To our knowledge, this is the first time that these two calculational 
approaches have been directly compared, therefore the agreement between them is noteworthy. In addition, we 
interpret this agreement as a very robust test of the veracity of our results. 




FIG. 1: Feynman diagram representing the pcrturbative contribution to the correlation functions. Solid and curly 
lines represent quark and gluon propagators respectively, while the dashed line represents the interpolating current. 




FIG. 2: Feynman diagram representing the {agGG) contribution to the correlation functions. 

The correlation function obtained at the quark-gluon level should be equivalent to that described at the hadron level 
due to quark-hadron duality, establishing a sum rule to extract the hadron mass. The Borel transform is applied to the 
correlation functions at both levels to pick out the lowest lying resonance, eliminate dispersion-relation subtractions, 
and enhance the OPE convergence. Using the spectral function in Eq. ([S|), we arrive at 



Ck{so,Ml)=f_ 



/x'^x 



S+2fc^-m^/A/| 



ds s ~ p{s) e 



-s/Ml 



(10) 



where sq is the continuum threshold parameter and Mb is the Borel mass. Then the hadron mass can be extracted 
using 



'X 



Ci (50 , Ml) 
Co {so, Ml) 



(11) 



in which mx denotes the heavy hybrid mass. 

For all the interpolating currents in Eq. ([2]), we calculate the correlation functions and the spectral densities are 
then obtained via p{s) = — Imn(s). We list the expressions for the correlation functions and the spectral densities in 
Appendix [A] As mentioned above, we have used two distinct methods to determine these expressions for the spectral 
functions. First, we have used the same method as in Refs. [3a, [Sa] to calculate the full correlation function, from 
which the imaginary part can be extracted through analytic continuation and hence a closed form expression for 
the spectral function can be constructed. Second, we have used the method of Refs. [60l - l63| which can be used to 
determine integral representations of the spectral functions. We have numerically verified that these two approaches 
lead to identical results. This provides a very strong test of our results. Another check is that our results for the 
perturbative and gluon condensate contributions are in numerical agreement with those in Refs. |30h32I |. 

Besides the perturbative term and the dimension four gluon condensate (asGG), we also calculate the dimension six 
condensates (g^fGGG) and (g^jj)- The condensate {gtij) has not previously been evaluated in the hybrid sum rules. 
As mentioned above, it comes from applying equations of motion in the condensate (DDG) arising from the diagram 
Fig. ((5^). Within the vacuum factorization assumption, (g^jj) is proportional to the square oi g^{qq). Although the 
numerical analysis shows that (g^jj) is small enough to be ignored compared to (g^/GGG), we include its contribution 
in our analysis because of the important stabilizing role of the dimension six condensates. We have only included the 
integral representations for these contributions in Appendix [^ 

From the formulae for the spectral densities in Appendix 13 we find that the perturbative contributions are in- 
variant and the gluon condensate and tri-gluon condensate contributions change sign under parity reversal. This is 




FIG. 3: Feynman diagrams representing the (DDG) (a) and 

functions. 



g^fGGG) (b) contributions to the correlation 



in agreement with the rcsuhs for the perturbative and gluon condensate contributions in Refs. |3Cll - l32| . It should be 
noted that both the expression for p^'^^^') and that of p^^-'^ have a singularity at s = 4m^. One should be very cautious 
in calcula ting the Borel transforms of the correlation functions from these expressions due to these singularities. In 
Refs. J35l |36| a numerical limiting procedure was developed based on the full analytic structure of the correlation 
functions and the relation between the Borel transform and inverse Laplace transform. This limiting procedure serves 
to cancel the integration divergences at s = 4to^ in the 77 ^ limit. A suitable value of rj should be determined to 
perform the numerical analysis. Alternatively, these divergences can also be eliminated if we use the integral forms of 
p(GGG) Q^^ pijj) jj^ gqg^ (|A8[) - (|A16I) where the delta functions absorb divergences of the spectral densities at a; = 
and X = 1. The 77 — > limiting expressions are necessary when the closed form expressions for the spectral densities 
are used (see Appendix IX)) . 

III. NUMERICAL ANALYSIS 

To perform the QCD sum rule numerical analysis, we use the following values of the heavy quark masses and the 
condensates [641 - 1671 : 

mc{^J. = TUc) = TOc = (1-28 ± 0.02) GeV, 

mb{n = rrib) = nib ^ (4.17 ± 0.02) GeV, 

(asGG) = (7.5 ± 2.0) x lO^^ GeV^, 

{glfGGG) ^ -(8.2 ± 1.0) GeV^(asGG), (12) 

{qq) = -(0.23 ±0.03)^ GeV^ 

in which the charm and bottom quark masses are the running masses in the MS scheme. Note that there is a minus 
sign implicitly included in the definition of the coupling constant g,, in this work. The scale dependence of these MS 
masses can be expressed in the leading order expressions: 

/ f \ \ 12/25 

■md^i) ^mA — T^TT- ) , (13) 



as{mc) 



12/23 



mfc(^) = mfcf^^) . (14) 

For the hybrid charmonium and bottomonium analyses, the strong coupling is determined by evolution from the r 



and Z masses, respectively: 



as{Mr) 



asip) = 25a (M ■) , — r;;^' ^--^i^^r) = 0.33; (15) 



as{Mz) 



aM = ,0^ I, C —, asiMz) = 0.118. (16) 

-^ ^ 127T ^"SVXFJ'' 

in which the r and Z masses, as{Mr) and as{Mz) are from the Particle Data Group 0. 

We begin with the sum-rule analysis of the hybrid charmonium mass spectrum. The stability of the QCD sum 
rule requires suitable working regions for the continuum threshold parameter sq and Borel mass Mg. According 
to Eq. (Uni), excited state contributions are naturally suppressed by the exponential weight function for small Mg. 
The lowest lying resonance contribution will be enhanced in the same region. On the other hand, however, the OPE 
convergence suffers if Mg is too small. In our analysis, the Borel mass working region is determined by the convergence 
of the OPE series and the pole contribution. The requirement of OPE convergence determines the lower bound on 
Mg while the pole contribution constraint leads to the upper bound. 

For hybrid charmonium the dominant nonperturbative contributions come from the gluon condensate (au GG) . 
However, the dimension six condensates are also important because they can stabilize the mass sum rules |34l436l | . 
We have calculated dimension six condensates (glfGGG) and {g^jj) in the spectral densities, extending the results 
of J30l - l32| . From the numerical analysis, we see that the condensate (^fjj) is much smaller than {glfGGG) (Fig.|l|). 
To study the OPE convergence, we require that the gluon condensate be less than one third of the perturbative term 
while the tri-gluon condensate be less than one third of the gluon condensate. This requirement gives the lower bound 
on M^. For example, we show the OPE convergence for the exotic channel J^*^ = 1 ^ in Fig. 21 from which we 
obtain the lower bound of the Borel mass Af^„„;„ = 4.6 GeV^. 
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FIG. 4: The contributions of each term in OPE series, including the perturbative term, gluon condensate {agGG), 
tri-gluon condensate (g^fGGG) and the condensate (^fjj) for the 1 ^ charmonium hybrid sum rule Co{Mg) with 
So —^ oo. 

We define the pole contribution (PC) as: 



PC(so,Af|) = 






(17) 



where PC is a function of sq and Mg. In order to study its dependence on A/^. one should first fix the value of sq. 
In order to do so we study the variation of the hybrid mass mx versus Sq, as shown in Fig. [5] for the J^^ = 1 ^ 
channel. We find that the hybrid mass is very stable for different values of M^ around sq = 17 GeV^. Using this 
value of So, we obtain the upper bound on the Borel mass of M^^^^ = 6.5 GeV-^ by requiring the pole contribution 
be larger than 10% in Eq. pT)) . In Fig. [5l we plot the Borel curve of tox versus Mg for the 1 ^ charmonium hybrid, 
which shows that the mass sum rule is very stable in the Mg working region. Finally, the predicted mass of the exotic 
J^'-' = 1 ^ hybrid charmonium is 3.70 GeV. This value is about 0.5 GeV lower than the lattice result in Ref. [y]. 
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FIG. 5: The variations of mx with so and Alg for the J^'^ = 1 + charmonium hybrid. 



After performing the QCD sum rule analyses, we extract masses of the charmonium hybrids with various quantum 
numbers as summarized in Table HI We also give the corresponding Borel windows, threshold values and pole con- 
tributions. Only errors from the uncertainties in the charm quark mass and the condensates are taken into account. 
We do not consider other possible error sources such as truncation of the OPE series, the uncertainty of the thresh- 
old value So and the variation of Borel mass Mb- In fact, the dominant error is from the uncertainty of the gluon 
condensate for the channels with J^^ = , 0+"'", 1"''"'', 1'' , 2++ while the tri-gluon condensate dominates for those 
with J^'~^ = '", 0"' , 1 '', 1 , 2 ^ . The errors from the charm quark mass rric and {g'^jj) are much smaller. 



J^^' 


so(GeV^) [M-^.^ML^^]{GeV') 


mx(GeV) PC(%) 


1-- 


15 


2.5- 


-4.8 


3.36 ±0.15 


18.3 


0-+ 


16 


5.6- 


-^7.0 


3.61 ±0.21 


15.4 


1-+ 


17 


4.6- 


-6.5 


3.70 ±0.21 


18.8 


2-+ 


18 


3.9- 


-7.2 


4.04 ±0.23 


26.0 


0+- 


20 


6.0- 


-7.4 


4.09 ±0.23 


15.5 


2++ 


23 


3.9- 


-7.5 


4.45 ±0.27 


21.5 


1+- 


24 


2.5- 


-8.4 


4.53 ±0.23 


33.2 


1++ 


30 


4.6- 


' 11.4 


5.06 ±0.44 


30.4 


0++ 


34 


5.6- 


' 14.6 


5.34 ±0.45 


36.3 


0— 


35 


6.0- 


' 12.3 


5.51 ±0.50 


31.0 



TABLE I: Masses of the charmoniuni hybrid states and the correspondmg sg, Borel windows and pole contributions. 

~7^' ,so(GeV^) [M^-^,A/fkax](GeV^) mx(GeV) PC(%T 



1-- 


105 


11- 


- 17 


9.70 ±0.12 


17.2 


0-+ 


104 


14- 


-16 


9.68 ±0.29 


17.3 


1-+ 


107 


13- 


-19 


9.79 ±0.22 


20.4 


2-+ 


105 


12- 


-19 


9.93 ±0.21 


21.7 


0+- 


114 


14- 


-19 


10.17 ±0.22 


17.6 


2++ 


120 


12- 


-20 


10.64 ±0.33 


19.7 


1+- 


123 


10- 


-21 


10.70 ±0.53 


28.5 


1++ 


134 


13- 


-27 


11.09 ±0.60 


27.7 


0++ 


137 


13- 


-31 


11.20 ±0.48 


30.0 


0~ 


142 


14- 


-25 


11.48 ±0.75 


24.1 



TABLE II: Masses of the bottomoniuni hybrid states and the corresponding sq, Borel windows and pole contributions. 



The results show that for the unstable channels with J^'-^ = '', 0^ , 1 ^, 1 , 2 ^ in Refs. p3C)| - l33 |. the dimension 
six condensate (g^fGGG) in the spectral functions can stabilize the systems and make it possible to extract the hybrid 
masses. This universal trend confirms the stabilizing effect of dimension six condensates previously observed for the 
J^*-^ = 1 [3J], 1++ [Sal and ^ [3g] channels. For these previously unstable channels, the pole contributions are 
very small even though we consider the dimension six condensates contributions. The small pole contributions result 
in the narrow Borel windows compared to the stable channels in Refs. [33432 with J^*^ = , O"*"*", 1^"'", 1^ , 2++. 
The 1++ and '' results agree with the previous sum-rule analyses including the dimension six gluon condensates 
[3a,|3g|. However, the 1 mass is somewhat smaller than in Ref. [SJl a result we attribute to the absence of the QCD 
continuum from the condensates in Ref. |34l |. We obtain four charmonium hybrid states with J^ ~ (0, 1, 2) ^, 1 
in the range 3.4 — 3.9 GeV which are much lower than the other channels in Table HI However, the 1 channel 
was examined in Ref. (68[, finding a mass of 3.66 GeV, which is comparable to our prediction for this channel. In 
the MIT bag model [j, Q , the light hybrids with these quantum numbers have been predicted to form the lightest 

PC _L 

hybrid supermultiplet consisting of a S-wave color-octet qq pair coupled to a Jg " " = 1"' gluonic excitation. This 
supermultiplet was confirmed in the heavy quark sector using lattice QCD [y| and the P-wave quasigluon approach j43 |. 
Our result supports this prediction although the extracted masses are 0.4 — 0.7 GeV lower. A heavier hybrid 
supermultiplet in Refs. [a, |43 contains states with J^'-^ = ^, (1 ^)^, (2 '')^,3 ^, (0, 1,2)++. It is composed of a 
P-wave color-octet qq pair coupled to a gluonic field with Jg ' " = 1^ . In Table U we list our mass predictions for 
the J^'^ = (0, 1)^ , (0, 1, 2)++ members of this excited hybrid supermultiplet. 

We obtain three charmonium hybrid states with exotic quantum numbers J^'~^' ~ 1 ^,0^ ,0^^ in Table HI The 
experimental identification of these states is considered to be a smoking gun for the existence of the gluonic degree of 
freedom in QCD because they cannot mix with the conventional charmonium states. The lightest exotic charmonium 
hybrid state with J — 1 ^ in the first hybrid supermultiplet is particularly important. The mass of the heaviest 
hybrid charmonium with J^'^ = 0^^ is about 5.5 GeV, which is consistent with the lattice QCD prediction [33| 
within the errors. This large mass may suggest that it has highly excited gluonic structures. 

In principle, conventional charmonium states can also couple to the non-exotic charmonium hybrid currents. Al- 
though a full analysis of the mixed hybrid/charmonium scenario is beyond the scope of this work (see e.g., Ref. |68| ) 
we can explore the qualitative effect of mixing on our mass predictions by modifying the right-hand side of (jlOp to 



include a conventional charmonium ground state rricc 

Cu (.0, Ml) = /im^+2'=e-"-/^^« + /^m^+^'^e"™-/^^: 

so that pT|) becomes 

2 _£i(,5o,M|)-/2miOe 



ds s'' p{s) e-"/^^« , 



l-JMl 



C,[s,,Ml)-P^ml,e-^l^J'^l 



(18) 



(19) 



By inputting the known charmoniuni mass nice and allowing fc to increase from zero (i.e. the pure hybrid case) we 
can determine how the hybrid mass is influenced by mixing. Fig. [B] shows that the 1 hybrid state mass increases 
as mixing with conventional charmonium increases and similar behaviour is found for all the other non-exotic cases. 
Mixing effects will thus tend to raise the mass predictions of Table |T] implying that our results are a lower bound on 
the mixed hybrid mass. A conservative estimate would be a mixed mass in the upper range of the uncertainties given 
in Table H 
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FIG. 6: The variations of mx with fc for J^'-^ = 1 hybrid charmonium system. The Borel scale is the central value 
from Table U and we have chosen several values of sg near the optimum value of Table HI 

The quantum numbers of X(3872) are given as J^'^ = 1++ or 2 ^ in the PDG Q- Although the analysis of 
angular distributions favors the assignment J^^ = 1++ [6^, [ZO], the 2 ^ assignment is also possible iTll. Recently, 
the LHCb collaboration has provided strong evidence for J^'-^ = 1++ UM- ^^ previously found in Ref. [Sa], the mass 
of the l^"*" hybrid charmonium in Table|T]is around 5.06 GeV, which is much higher than the mass of X(3872), which 
seems to preclude a pure hybrid interpretation as pointed out in Ref. [35|- The mass of the 2 '" hybrid charmonium 
is about 4.45 GeV, which also precludes the 2 ^ hybrid charmonium assignment for the X(3872), in agreement with 
the LHCb result 
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hybrid bottomonium system. 



By replacing rric with mt in the spectral densities and using the strong coupling in Eq. (|16p . we perform the same 
analysis as described above and collect the numerical results for bottomonium hybrid states in Table |IT1 Obviously the 
Borel windows are enlarged compared to the hybrid charmonium systems, which means that the stabilities are better 
for the bottomonium hybrid states. We show the Borel curves of the bottomonium hybrid state with J-^'-^ = 1 ^ in 
Fig. [71 In Table HIl the masses of the four bottomonium hybrid states with J^*^ = (0, 1, 2) ^, 1 are about 9.7 ~ 9.9 
GeV. The mass variations of these states are much smaller than those in the charmonium system. They form the 
lightest bottomonium supermultiplet, in complete analogy with the charmonium hybrid supermultiplet structure. In 
Table HIl the masses of ^ and 1++ bottomonium hybrids are again in agreement with Refs. jS^ |36i] and the 
state is still the heaviest one. 



IV. SUMMARY 

In this paper we have studied the mass spectrum of charmonium and bottomonium hybrids using QCD sum rules. 
Using the hybrid operators in Eq. ^, we calculate the correlation functions and the spectral densities, confirming the 
perturbative and gluon condensate results of |30l432l | and extending them to include the dimension six condensates 
{glfGGG) and {g^jj). Agreement with the (g^fGGG) results of Refs. [Sl, |3g| provides an additional consistency 
check on our analysis. Expressions for the correlation functions and spectral densities are given in Appendix \X\ 

The Feynman diagram in Fig. ([3^) gives contributions to both the condensates (g^/GGG) and {g'^jj). The conden- 
sate (g^jj) has not previously been calculated in hybrid sum rules. Given the crucial stabilizing contribution of the 
dimension six gluon condensate, we have also included the dimension six quark condensate in our analysis. The chan- 
nels with J^*-^ = ^, 0"' , 1 '", 1 , 2 '" were unstable without the contributions of the dimension six condensates 
in the original QCD sum rule analysis J3(]| - l32l |. However, these channels are stabilized by the effects of the dimension 
six condensates in our analysis. In our analysis all channels lead to stable mass sum rules. Due to the small values 
of the pole contribution, the Borel windows of the channels with J^^ = ^, 0"' , 1 '", 1 , 2 ^ are much narrower 
than those of J^^ = ,0++,l++,l^ ,2++. We speculate that this explains why these channels were unstable 
in Refs. |30l432| : the pole contributions were so small that no effective working regions for Borel parameter existed. 
Although the {g'^jj) effects are very small compared to (g^fGGG), we still include its contribution to the correlation 
function in this work because of the important stabilizing role of the dimension six condensates. 

The sum-rule mass predictions show that for the charmonium hybrids, there are four negative-parity states with 
J^'-^ = (0,1,2) '",1 which lie below 4 GeV. Their masses are about 3.4 ~ 3.9 GeV and are much lower than 
the other states, forming the lightest charmonium hybrid supermultiplet. Our mass prediction for the 1 state 
is comparable to the result found in Ref. [6a], which provides indirect support for our mass predictions for the 
states in the lightest supermultiplet. These states were predicted to form the lightest hybrid supermultiplet in lattice 
QCD [6[ and the P-wave quasigluon approach [42[, consisting of a S-wave color-octet qq pair and a 1"' excited gluonic 
field. A heavier hybrid supermultiplet is composed of a P-wave color-octet qq pair coupled to a 1"' excited gluonic 
field [y, \^ . In our spectrum, we predict masses of five positive-parity members of this excited hybrid supermultiplet 
with J^*-^ = (0, 1)"' , (0, 1, 2)++. We speculate that the supermultiplets are occupied by states of the same parity 
because of the behaviour of the perturbative, gluon condensate and dimension six condensate contributions under 
parity reversal. In particular. Appendix |^ shows that the perturbative contributions are invariant and the gluon 
condensate and tri-gluon condensate contributions change sign under parity reversal, implying a consistent pattern in 
the mass hierarchy of the positive- and negative-parity states. We note that the J^'-' = hybrid has a much higher 
mass in our spectrum which may suggest that it involves a different excitation of the gluonic field compare to other 
channels. We list the masses of the charmonium and bottomonium hybrid states in Table Hand Table HIl respectively. 
Finally, we have investigated the effects of mixing between non-exotic hybrids and conventional quarkonia. We find 
that mixing increases the mass predictions for the non-exotic hybrids given in Tables U and |TI1 and we conservatively 
estimate that the mixed states masses are in the upper range of the uncertainties given in Tables U and |TT1 A more 
detailed study of mixing between non-exotic hybrids and conventional quarkonia is left for future work. 

To date, there are many interpretations of the X(3872) meson, such as a molecular state, a tetraquark state, and 
a charmonium hybrid. LHCb has provided strong evidence for the J^^ — l^'^ quantum numbers, seemingly ruling 
out the 2 ^ assignment [T^. In our spectrum, the masses of the 1++ and 2 '" charmonium hybrids are about 5.06 
GeV and 4.45 GeV, respectively. The 1"'"+ (also see Ref. [33) and 2 ^ mass predictions seem to preclude a pure 
charmonium hybrid interpretation of the A'(3872) meson. 

Observation of the overpopulation of the states with qq quantum numbers in the charmonium and bottomonium 
regions, such as the X, Y, Z states, is an important signal for the existence of the heavy quark hybrids. As an 
experimental signal, this overpopulation complements the identification of states with exotic quantum numbers. 
Experiments such as BESIII, PANDA and LHCb will collect more definitive data on the charmonium spectrum, 
including hybrid mesons. Our results for the mass spectra of charmonium and bottomonium hybrids in Table [J 
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and Table |ll] may be useful for the future search for these fascinating exotic and non-exotic hybrid states at these 
experimental facilities. 
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Appendix A: Correlation Functions and Spectral Densities 

In this appendix, we list the correlation functions and spectral densities of the hybrid interpolating currents in 
Eq. ^. First, we tabulate the correlation functions obtained by the first approach mentioned in Sec|lTl For all of the 
various J^^ quantum numbers studied in this paper, the perturbative contributions can be expressed as 
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(^) 



mPas 
16207r3 
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i,i,i;2'3;^ 



f2iz)3F2 



l,l,2;-,4;z 



./aW 



4j^2 



(Al) 



where we have omitted terms corresponding to dispersion relation subtraction constants, 3F2 denotes the generalized 
hypergcomctric function [73| and the functions /i (z), /2 (z) and f^ (z) arc polynomials in z that are tabulated in 
Table Hill for each heavy hybrid channel studied in this paper. For the dimension four gluon condensate (agGG), all 
contributions can be expressed as 
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where we have again omitted terms corresponding to dispersion relation subtraction constants, 2-F1 denotes the Gauss 
hypergcomctric function [73l |. z is defined as in (jAip and the function /4 [z) arc polynomials in z that arc tabulated 
in Tabic Hill for each heavy hybrid channel studied in this paper. 
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TABLE III: Polynomials /i (2), /2 (2), /s (2) for the perturbative contribution ()A1|) and /4 (2) for the gluon condensate 
contribution (|A2p . Note that for J^'~" = 1 ^ , 0"'""'" , 1"' , , 2++ omitted above, the polynomials /i (2), /2 (2), f^ (2) 
are identical to those for J^^ = 1++ , ^ , 1 , 0"' , 2 ^, whereas the polynomials /4 (2) differ by an overall minus 
sign. 



Now we consider contributions from the dimension six gluon condensate, (g^GGG). For convenience, we tabulate 
contributions from the left diagram and right diagram in Fig. [3] seperately. The left diagram, which arises from the 
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vacuum expectation value (DDG) is denoted as Hqqq, while the right diagram is denoted as Hqqq. All contributions 
to IIqqq can be expressed as 



n£oc(^)= ^1(f,;f\,_\.. (/5(^) + /a(2).F, 



1'1'2=" 



Similarly, all contributions to Hqqq can be expressed as 
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where in Eqns. (jA3p and (|A4p we have omitted terms corresponding to dispersion relation subtraction constants, all 
notations are identical to those in (|A2p and the functions f^ (z), fe {z), /? (z) and /§ (z) arc polynomials in z that are 
tabulated in Table ITVl for each heavy hybrid channel studied in this paper. 



J^^ 


hiz) 






hiz) 




f7iz) 


/8(^) 


1++ 


51- 138z + 81z^ 


- 


22(2 


-92 


+ 622) 


27-512 


42(-l + 32) 


0-+ 


-48 + 1502-93^2 


-3 + 62 


-182 


2 + 1223 


-27 + 632 


-1222 


1-- 


78- 1562 + 812^ 


-3 


- 42 + I82 


2 - 1223 


63-512 


42(-4 + 32) 


0+- 


-3(-l + 2)(-25 + 3l2) 


6(- 


-1 + 


z)z{- 


1 + 22) 


63(-l + 2) 


-12(-1 + 2)2 


2-+ 


i (5 - 462 + 322^) 






1 

2 










1-+ 


-3 (-3 -62 + 722) 










3(3 + 52) 




0++ 


-12-302 + 332^ 










-9 - 272 




1+- 


-3(6-122 + 722) 










-3(9-52) 




0— 


-3(5-ll2)(-l + 2) 










-27(-l + 2) 




2++ 


i(-13 + 622-4022) 
















TABLE IV: Polynomials /s (2), f^ (2) for the contribution (jA3p and fr (2), /§ (2) for the contribution (jA4p to the 
dimension six gluon condensate. Note that the polynomials fe (2) and /s (2) for J^*^ = 1 ^ ,0++ , 1"' ,0 ,2++ 
only differ from those for J^'^ = 1++ , ^ , 1 , 0^ , 2 ^ by an overall minus sign. 

We do not include explicit expressions for the dimension six quark condensate correlation function, although integral 
representations for these contributions are included in the following. 

Note also that the imaginary parts corresponding to the perturbative (|Al[) . dimension four gluon condensate (jA2p 
and dimension six gluon condensate (|A3|) , (jA4p correlation functions can be extracted through analytic continuations 
of the hypergeometric functions in these expressions. This is precisely what was done in Refs [3^, l3g|- Given the 
universal forms of these hypergeometric functions for all J^'~' channels, we could construct tables similar to those 
above for each imaginary part. However, we have not chosen to do so. Instead, we have listed explicit expressions for 
these imaginary parts in the form of spectral densities in order to emphasize the equivalence of the two approaches 
that we have taken to calculate the Wilson coefficients. In what follows, all spectral densities will be given in terms 
of a closed form expression and an integral representation, and we have verified that these are numerically identical. 

Using the second approach mentioned in Sec. |TT1 we calculate the spectral densities up to dimension six: 

p(s) = /-* (s) + p<GG) (^) ^ piGGG) (g) ^ pin) (g) ^ ( A5) 

The integration limits in the following expressions are: 
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y 1 — Am? / s^ 



CoiM^B^so 



So 



/■So 

J4m2(l+T;) 






e-'*/^-fSrfs+ lini 

T;-i.O+ 



■ ds 



4m^ 

7^ 



1287r2 



967r2 



(AS) 



• For J^^ = 1 



/^'■*(s) = 


pr-+(^)> 










p<««)(.) = 


-pf-f{sl 




p<«««>(s) = 


-pf-T\^l 




^(«>(s) = 


-2pi^ii(s), 




Co{M%,so) ^ 


J4m2 


/^-(.)+p<««>(.y 


e"'/^^^(is+ lim \ 






JAm^il+ri) 


p<'^'^'^>(s)+pto>(s)' 


-^ 4m2 


lig^GGG) {gim 
1927r2 1447r2 


. For J^C7^o-+. 


















p<^'^> (s) 



{gin) f\ r(3-8x2 + 4x3)m4 . 







2887r2 



-2x0{s-m^)\ 






2(2.t3 _ 5a;2 + 6a; - 4)m^ 



(A9) 



c(l-a;) 



(l-x)3 



(5(s — to2) 



(AlO) 



487r2(s-4m2)2s 



a/I - 4m2/s, 



Co{Ml,so) 



L 



So 

So 






e-^A4ds+ lim 



Am"^ {1+r}) 



e "Bds + 



4m2 



{glJGGG) , (g^jj) 



1287r2 1927r2 



For JP^ = 1+- 



pP'='"*(s) 



da 



Pma. 



dp 



^2(1 - a - /3) f (a + ;3)3(9a2 + 12a;3 - 9a + 3/32 _ ^^^ 



967r4 V 



2«3 



a2/3' 



3(a + /3)2(20a2 + 24a/3 - 25a + 4/32 „ g^ + 5)^4^ 
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+ 



3(a + /3)(35a2 + 40a/3 - 49a + 5/^2 - 21/3 + 14)to2s2 

a{bAa^ + 60a;3 - 81a + 6/?^ - 37/^ + 27)s^\, 

15 log [^/i + \/7^ (-1 + 12z - 48z2 + 1282^) 



7207r3z2 



V^V^^l (-15 + 170z - 608^2 - 1104z3 + 192z'' 



367r 

1927r"^ Jo '^ 2;(1 — x')'^ 
-(l + 2a;)e'(s-m2)| 



4to2 



3-1^2 



(1 -8a; + 4a;2 + 2a;'^)m' 
a;(l-x)2 



■5{s — m ) 



m^{glfGGG) hs^ - 24m2s + 24™" 
967r2 (s — Arri^Ys 



■\/l — 4to2/s, 
(1 - Ax)m? 



(gtJJ) 3s^ - ISm^s^ + 8m^s + 96m^ r- 

288^ (s - 4m2)2s ^1 - 4m /s, 



/:o(A/|,so) 



So 

so 









4m2(l+rj) 



e ''^B ds 



7^ 



(g3/GGG) , (g^jj) 



3847r2 2887r2 



(All) 



• For J^^ = 



.pert/ \ _ 



pP"'''(s) 



rfa 



/3„„x g2(^ _ a - /3) r (a + /3)3(9a2 + 12a/3 - 9a + 3/32 _ ^^^6 



d^- 



/3„i„ 



967r4 



■{ 



2«3 



a2/3' 



3(a + /3)(12a3 + 36a2^ - 15a2 + 3a - 24a/3 + 36a;32 _,_ j^2/?3 _ 9^2 _ ^^^ 
3(15a3 + 55a2/3 - 21a2 + 6a - 46a/3 + 65a/32 + 25/^3 - 25/32 + 2l3)m^s^ 

-a(18a2 + 60a/3 - 27a + 42/3^ - 49/3 + 9)s3|, 

fi / 

15 log [yi + V^^ (-1 + 6z - 16^2 _^ 32^3^ 



2407r3z2 



-^/I\/I^T (-15 + 80z - 188^2 - 224z^ + 322"*) , 



4m2 



,<««)(.) = - <"-^^^^(;^-'"^'Vl-4m2/. , 



JGGG) 



(5) 



(g3/GGG) /•! r2(l-2x) 



^r-{^ 



1927r2 7o " L a;(l-a;)3 
-{l + 2x)e{s~m^)\ 
n?{glfGGG) 



m'^ ^,, _2n 2(1 - 5x + 2a;2)m2 ^„ 

S {s - m^) ^ -^ 6{s - m^) 



x{l — x) 



327r2(.s-4m2) 



■\/l — 4to2/s, 
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pto)(,)^i^ /' dx{ ^"-^ ~ 'y^ -^Z ~ '''" 3'i^ - ^') 



Ix^ - 16x2 + 8a;- l)m'' 



2887r2 Jo I x{l ~ x)4 

2(8x3 - 22x2 ^ 22x - 7)m2 

{g^ij) s^ — 2m2s + 8m^ 



(5(s-m2) + 4x6i(s-7n2)| 



(A12) 



Co{Ml,so) 



967r2 (s — 4r7i2)s 

So r 



■\/l — Am? / s^ 



e "b ds + lim 



p^^^^^(s)+p^-'^'(s) 



O'i) I 



e *^B ds, 



• For JP^ = l—-. 

/-*(s) = p^:i*(s), 



pi^^^Hs)^-pf,^^\s), 



{gtij) f^ , r(3-16x + 16x2)m4 ^^ (2x2 - 4x + 5) 



""•'^''-iHiil-^^ .(I-.). 



S'{s ~ m2) + 



(l-xy 



> 

'-5{s^m^)\ (A13) 



tjj) 3s3 - 18m2s2 + 20m'^s + 24m^ 



{91J3 



A (Ml, so) 



1447r2 

So 



(s - 4m2)2 



^1 - 4to2/s, 



4m2 
So 



4m2(l+r;) 



pP^'-*(s)+p<'=«>(s) 

<GGG)(s) + pto>(s) 



e"''/^^Bds+ lim 



e *^B ds 



4m2 



2 r/„3 



{g^GGG) , (5^jj) 



3847r2 5767r2 



For JP^ = 0^ 



p^-*(s)=pr-(^)> 
p'^^Ks)^-pf_^\s), 



pin)(^s) 



{gtij) f\, r (16x3 - 20x2 + 4x + l)m4 



dx 



x(l-x)'' 



5 [s — rh ) 



2(2x3 - 7x2 _,_ j^Q^, _ 4)^ 

(g.'jj)s'-2m2s-4m4 , 

487r2 (s-4m2)s ^' *"'/''' 



-(5(s-m2) + 2x6l(s-m2)| 



/■So p 

£o(M|,so)=/ p^^'-*(s)+p<^^>(s) 



e "b ds + lim 



(A14) 



p^^^^^(s)+p^-'^'(s) 



(ii> ^ 



i ds. 



• For J-f"^ = 2-+: 
pP'='^*(s) 



2 03 



a2/3 



""^ , /■''""'== ,„fl2M _ Q, _ ^) , (^ _^ /3)3(6a2 + Tap - 6a + /?2 - S/?)/?!*^ 
Q{a + ;3)2(8a2 - IQq + 9q;/3 + ^2 _ 5^ _^ 2)m^s 
9(a + ^)(10a2 - 14a + lla;3 - 7/? + /?2 + 4)to2s2 



a/3 



-4a(12a2 - 18a + 13a/^ - 9/3 + /32 + 6)s3\ 
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m^as 



345607r3z 



3,3 



15 log [y/I + VI^T\ (7 - 256z^ + 11522 



+yi\/I^T (105 + 702 + 56^2 - 3792z^ - 125442"* + 25602^) 

^ ^ ^" 96^2 y^ °'''(l-x)2''^* "" ' 

-yl — Aim? / s^ 

,(..)(,). _^^/ dx{^^ 



(g3/GG'G)m6 
247r2(s - Am^Ys 

{gjjj) [' 

5767r2 io 



da;<| — -t: f;^ — 5'{s - m^) + ^ -—^ 5(s - m^) 



(1 - xf 
■{Ax^ - 2x)e{s - m^)\ 

{gtjj) ■s^ ~ Gm'^s'^ + 6m^s + 32to^ /- — ^i" 

'864;? is - 4m2)2s VI -4m /s, 



(l-x)2 



>Co(M|,so) 



So 
So 



4m^{l+Ti) 






e~^/*^«ds+ lim 

7;->0+ 



4m2 ' 



e "sds- 



Am'^ 

7^ 



(g^/GGG) , (.gfjj) 



7687r2 



11527r2 



e -B 



(A15) 



• For JP^ = 2^ 






(5;*JJ) /•'.r 2(3 -2a;) 



/fill 
I fl 



^'^'(^)=-|^X ^H ":rtm ^-^^)+ 



(4a;2 - 6x + l)m^ 
(l-x)2 



J(s — 771 ) 



-(4a;^-2a;)6'(s-7n^)| 



A (Ml, so) 



8647r2 

So 



4m2 
So 



(s — 4?Tt2)2s 

pP"=''*(s) + p<^^>(s) 



Vl - 4m2/i 



e-"/^^«ds+ lim 






ds 



4771"^ 

7^ 



(.gf/GGG) , {gtn) 



7687r2 



5767r2 



(A16) 



in which m is the heavy quark mass, 5' {s — rh ) 



2\ _ dS(s-m~) ~ 2 _ m 



as ' 



!"_ . where a; is a Feynman parameter. 



a:{l — a:) 



